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ABSTRACT 

We report the discovery of a gravitationally lensed hyperluminous infrared galaxy (intrin¬ 
sic I/iR ^ with strong radio emission (intrinsic I/i.4ghz ~ lO^^WHz”^) at 

z = 2.553. The source was identified in the citizen science project Space Warps through 
the visual inspection of tens of thousands of iJK^ colour composite images of Luminous Red 
Galaxies (LRGs), groups and clusters of galaxies and quasars. Appearing as a partial Einstein 
ring (re ~ 3") around an LRG at z = 0.2, the galaxy is extremely bright in the sub-millimetre 
for a cosmological source, with the thermal dust emission approaching 1 Jy at peak. The red- 
shift of the lensed galaxy is determined through the detection of the CO(3 -a 2) molecular 
emission line with the Large Millimetre Telescope’s Redshift Search Receiver and through 
[Olll] and Ha line detections in the near-infrared from Subaru/IRCS. We have resolved the 
radio emission with high resolution (300-400 mas) eMERLIN L-band and VLA C-band imag¬ 
ing. These observations are used in combination with the near-infrared imaging to construct 
a lens model, which indicates a lensing magnification of ft ^ 10. The source reconstruction 
© 0000 RAS appears to support a radio morphology comprised of a compact (< 250 pc) core and more 

extended component, perhaps indicative of an active nucleus and jet or lobe. 
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1 INTRODUCTION 

Strong gravitational lenses are striking astrophysical laboratories 
that allow us to study the structure and distribution of matter in and 
around massive galaxies and clusters (Koopmans et al. 2006, Ruff 
et al. 2011, Sonnenfeld et al. 2013, Richard et al. 2014), cosmogra¬ 
phy through quasar time delays (Refsdal 1964, Tewes et al. 2013, 
Suyu et al. 2013), perform highly detailed investigations of early 
galaxies (e.g. Kneib et al. 2005, Stark et al. 2008, Swinbank et al. 
2009) and even detect the first galaxies (Kneib et al. 2004, Zheng 
et al. 2012, Zitrin et al. 2014). The number of confirmed gravita¬ 
tional lenses is still relatively low: only some 500 are confirmed 
to date. With diverse astrophysical and cosmological applications, 
surveys for new strong gravitationally lensed systems remain a use¬ 
ful endeavour. In this article, we present new results from a novel 
search for new gravitational lens systems using wide-field near- 
infrared and optical data. Lens searches in the visible bands are 
typically optimised for finding blue arcs, representing the images 
of star-forming galaxies in the early Universe (which are abundant). 
Searches at longer (near-infrared and beyond) wavelengths have the 
advantage of selecting dust-obscured or very high-redshift sources 
that could easily be missed in optical searches (e.g. Negrello et al. 
2010; Vieira et al. 2013). 

The VISTA-CFHT Stripe 82 survey (VICS82, Geach et al. 
in preparation) is a wide area near-infrared survey of the Sloan 
Digital Sky Survey (SDSS) Stripe 82 with CFHTAVIRCam and 
VISTA/VIRCAM, covering 130 square degrees to a 5a depth of 
approximately 22nd AB magnitude in the J and bands. Stripe 
82 is also covered by deep (24th magnitude) /-band imaging from 
the Canada-France-Hawaii Telescope (CFHT) in the CS82 survey 
(Moraes et al. 2014, Erben et al. in prep) with comparable image 
quality (sub-arcsecond), as well as the deeper coverage from SDSS 
itself across the equatorial stripe. The combination of moderately 
deep optical-near-infrared imaging over a wide area is ideal for 
searches for strong, ‘red’ gravitational arcs. After serendipitously 
finding a handful of promising lens candidates through visual in¬ 
spection, we decided to pass the dataset to Space Warps (Marshall 
et al. 2015, More et al. 2015). Space Warps is a citizen science 
project to discover gravitational lenses run as part of the Zooni- 
verse family^ of online projects (Marshall et al. 2015, More et al. 
2015). 

Gravitational lens discovery remains a labour-intensive pro¬ 
cess, even automated algorithms produce highly impure samples 
that require visual inspection for verification (Marshall et al. 2009, 
More et al. 2012, Gavazzi et al. 2014, Chan et al. 2014, Brault & 
Gavazzi 2014). Space Warps was conceived to mitigate this issue 
for gravitational lens searches in the era of large surveys which have 
the potential to increase current sample sizes by several orders of 
magnitude. Through an online interface, users (the general public) 
are presented with tens of thousands of images, each of which are 
inspected for evidence of strong gravitational lensing. A sophisti¬ 
cated training scheme using images containing realistic simulated 
lenses, built into the inspection interface, is used to simultaneously 
train users and calibrate the classification process (Marshall et al. 
2015, More et al. 2015). Humans are highly proficient at detect¬ 
ing the - often very subtle - features associated with gravitational 
lenses, and so harnessing the power of a natural neural network of 
several tens of thousands of ‘nodes’ is a particularly efficient ap¬ 
proach to this problem. 

The VICS82+CS82 search with Space Warps (actually the 
^ http://zooniverse.org 



Figure 1. The colour composite image of 9io9 that was presented to online 
visitors to Space Warps . Users were presented with a (random) sequence 
of 40" cut-out images where the blue channel represents z-band data from 
the CHFT Stripe 82 (CS82) survey (Erben et al. in prep) and the green and 
red channels are J and Ag-band imaging from the VISTA-CFHT Stripe 82 
(Vies82) survey (Geach et al. in prep). Images were centred on Luminous 
Red Galaxies (LRGs), groups and clusters and quasars. Volunteers were 
asked to identify potential lensed features, having been trained using simu¬ 
lated images. In this case, a red ring of radius 3" can be seen around the 
central galaxy (an LRG at 2 ; = 0.2). This system was independently identi¬ 
fied by several users and quickly became a high-confidence lens candidate. 

second incarnation of the project, the first was an optical search 
within the CFHT-Legacy Survey) was accelerated by being 
launched on Stargazing Live!, broadcast on 7th January 2014 to 
approximately 2 million viewers in the United Kingdom by the 
British Broadcasting Corporation. The goal of Stargazing Live! 
is to engage the public with astronomy-related activities ranging 
from amateur pursuits through to the communication of the latest 
scientific research, major missions and projects. The wide reach 
of this programme (which was run over three consecutive nights) 
makes it an exceptionally successful and efficient public engage¬ 
ment tool for promoting astronomy, and indeed science in general. 
Space Warps was phenomenally successful: by the end of the final 
programme on 9th January, over 7.5 million classifications of im¬ 
ages had been logged, such that each //As-band image in the sam¬ 
ple was independently viewed approximately two hundred times. 
Out of several tens of excellent lens candidates identified by the 
collective efforts of all participants of Space Warps, the best can¬ 
didate we have appears to be an active, dust-obscured source at 
2: 2.5. Dubbed ‘9io9’ (a shortening of its Space Warps identifier 

ASW0009io9^), the target appears as an incomplete red Einstein 
ring of radius approximately 3" around a Luminous Red Galaxy 
(LRG) at 2 ; = 0.2 (Eigure 1). 

In this paper we present follow-up observations and lens mod¬ 
elling of 9io9 that reveal the galaxy to be an intrinsically radio- 
and sub-millimetre-bright hyperluminous infrared galaxy (HLIRG, 
Lir > lO^^L©) at 2 ; = 2.553. In Section 2 we describe archival 

^ http://talk.spacewarps.Org/#/subjects/ASW0009io9 
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Figure 2. Images of 9io9 (aligned north up, east left): left-to-right: (a) Kg-band image. The small point source (G2) to the north of the central (leasing) galaxy 
G1 is not a background source, but likely a satellite of the leasing galaxy at 2 ; = 0.202 (b) VLA FIRST 1.4 GHz image showing the radio detection, peaking 
slightly SW of the central LRG, (c) higher resolution 1.4 GHz VLA imaging from the survey of Hodge et al. (2011), showing that the radio emission resolves 
into a ring conhguration that traces the near-infrared light, but with a clear peak at the SW, (d) eMERLIN 1.4 GHz map (300 mas beam) resolving two main 
radio components: an arc (A) and point source (B), (e) VLA 5 GHz image at similar resolution, again resolving the radio morphology. Our lens model (§3.2) 
indicates that B is the counterimage of the bright, compact component at the southern extent of A. 


and new observations of 9io9, in Section 3 we constrain the inte¬ 
grated spectral energy distribution of the galaxy and model the lens 
itself, and in Section 4 we conclude with a brief discussion of the 
likely nature of 9io9. Throughout we assume a fiducial cosmology 
where Qm = 0.3, Qa = 0.7 and Hq — 70kms“^ Mpc“^. 

2 OBSERVATIONS 

2.1 Radio 

Immediately upon the discovery of the lens candidate dur¬ 
ing Space Warps we cross-checked the position of the source 
(02‘'09“41.3", +00° 15'58.6") in other imaging surveys. The Very 
Large Array (VLA) FIRST 1.4 GHz survey (Becker et al. 1994) 
reveals a 4 mJy radio source at this position, but the angular reso¬ 
lution of FIRST is too coarse to determine whether the radio emis¬ 
sion originates from the background galaxy, or the lensing galaxy. 
Hodge et al. (2011) present higher resolution VLA (A and B config¬ 
uration, 9 — 1.8") 1.4 GHz imaging of Stripe 82 that partially re¬ 
solves the radio fiux into a ring that follows the .^fg-band light (Fig¬ 
ure 2). This motivated us to obtain higher resolution radio imaging 
of 9io9, which we describe here. 


2.1.1 eMERLIN L-band 

We obtained eMERLIN observations over 10.25 hrs between 1415 
UT 9th January 2014 to 0030 UT 10th January 2014 (during the 
television broadcast). The target source was interspersed with the 
phase reference calibrator (J0208-0047) with a cycle time of 10 
minutes with 7 minutes on the target and 3 minutes on the phase 
reference. 30 seconds of data were flagged at the start of each scan 
to remove corrupted data logged during telescope drives with a re¬ 
sulting total time on the target source of 6.72 hours. The data were 
divided into 8 intermediate frequencies (IPs) with a total observing 
bandwidth of 512 MHz at a centre frequency of 1.5185 GHz. Data 
calibration was performed from calibration scans on the resolved 
flux density calibrator 3C286 (1331-1-305) and the point source cal¬ 
ibrator OQ208 (1407-1-284). The flux density scale assumed for 
3C286 incorporated the latest coefficients from Perley & Butler 
(2013) and a standard model of spatial resolution of the source as 
observed by eMERLIN. These were used to derive flux densities 
for the point source OQ208 which has a rising spectrum across the 
observed band (IPl: 0.815 Jy, IP8: 1.208 Jy), with a value at the 


central frequency of 1.0078 Jy. Estimated errors in the flux den¬ 
sity scale are at the level of 5-10%. Residual delay corrections for 
the data throughout the run were derived from all three calibration 
sources. After performing phase and minor gain corrections derived 
from the phase calibrator, the target was imaged and several rounds 
of phase-only self-calibration applied. The data calibration and im¬ 
age processing was performed in AlPS. The final beam is nearly cir¬ 
cular with PWHM 300 mas, and the map r.m.s. is 16 jiiy beam“^. 


2.1.2 Karl G. Jansky Very Large Array C-band 

Data were obtained during 1st March 2014 using the 3-bit samplers 
and an advanced mode of the immensely powerful WIDAR corre¬ 
lator at the heart of the Karl G. Jansky Very Large Array (VLA, 
proposal 13B-458), with eight sub-bands in each of the Al/Cl, 
A2/C2, Bl/Dl and B2/D2 basebands. We used the VLA’s C-band 
receivers with the correlator instructed to record a band compris¬ 
ing 1024 X 62.5 kHz full-polarisation channels centred on the 
22.23508 MHz water line, tracking a barycentric redshift of 2.5529. 
We also recorded 31 wider bands to sample continuum emission, 
each comprising 128 x 1 MHz dual-polarisation channels, covering 
4.05-7.96 GHz contiguously, with a few small regions of overlap. 
To mitigate time smearing, since we used the widest (A) configu¬ 
ration of the VLA, the integration time was set to 2 seconds, which 
resulted in a total data rate of 62Gbhr“^, with 26 x 287 seconds 
(roughly 2 hours) spent integrating on 9io9, 27 x 31 seconds on 
the gain calibrator, J0215—0222, and 226 seconds on 0137-1-331 
(3C 48) to set the flux density scale. The data were calibrated using 
version 4.1 of CASA, starting with a pass through the initial stages 
of the VLA data pipeline, then flagging extensively by hand, then 
re-starting the pipeline where we had previously paused. CASA was 
also used for imaging, which required us to map and clean a num¬ 
ber of bright, nearby sources to reduce contamination of the region 
around 9io9 by strong sidelobes. The resulting image of 9io9, cen¬ 
tred at 5 GHz, has a 500 x 400 mas beam (PA = 7°), with an r.m.s. 
noise level of 3-4/iJy beam“^. 


2.2 Sub-millimetre and millimetre 

A large area of Stripe 82 has been mapped at 250/im, 350/im 
and 500/im with Herschel, as part of the The Herschel Stripe 82 
Survey (HerS, Viero et al. 2014). At the position of 9io9 there is 
bright sub-millimetre detection, with 5'25o = 826 zb 7 mJy, S350 = 
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912 zb 7mJy, S 500 = 717 zb 8mJy. The low resolution of Her- 
schel cannot spatially resolve the emission, but the sub-millimetre 
colours are consistent with the emission originating from the back¬ 
ground source at z ~ 2-3 (the lensing galaxy is not expected to be 
a significant source of sub-millimetre emission). In order to better 
constrain the shape of the far-infrared spectral energy distribution 
and determine the redshift of the lensed source, we obtained fur¬ 
ther sub-millimetre imaging and spectroscopy that we describe in 
the following. 

2.2.7 James Clerk Maxwell Telescope SCUBA-2 

We obtained a 30 minute observation of 9io9 with the SCUBA-2 
camera on the James Clerk Maxwell Telescope (JCMT) on 15th 
January 2014 as part of Director’s Discretionary Time (project 
M13BD12). SCUBA-2 records both 450/im and 850/im data si¬ 
multaneously. We employed the DAISY scanning mode, appro¬ 
priate for point sources, and the zenith opacity at 225 GHz was 
T 225 GHZ = 0.07-0.1 over the course of the observation. The data 
were reduced with the dynamic iterative map maker of the SMURF 
package (Jenness et al. 2011, Chapin et al. 2013) and optimally fil¬ 
tered for detection of point sources. We applied flux conversion 
factors (FCFs) of 491 JypW“^ and 537JypW“^ at 450/im and 
850/im respectively, which are the canonical FCFs derived from 
observations of hundreds of standard sub-millimetre calibrators ob¬ 
served by SCUBA-2 since operations began (Dempsey et al. 2013). 
An additional correction of 10% is made to the FCFs to account for 
flux lost during the filtering steps. The absolute flux calibration at 
850/im is accurate to within 15%. 

9io9 is detected with = 167 zb 4mJy (the JCMT beam 
at 850/im is 15'', and so this can be considered an integrated mea¬ 
surement). At the same position in the 450/im map, we measure 
*S '450 = 217 zb 72 mJy, this is a factor of three lower than the Her- 
schel measured flux at approximately the same wavelength. Given 
the ~5x finer resolution of JCMT/SCUBA-2 compared to Her- 
schellSFlRE at 500/im, it is possible that SCUBA-2 is resolving the 
emission into several beams, rendering some of the sub-millimetre 
emission undetectable. This offers a first insight into the distribu¬ 
tion of cold dust in 9io9: higher resolution sub-millimetre mapping 
with a powerful interferometer such as the Atacama Large Millime¬ 
ter Array (ALMA) will allow us to compare the dust and radio mor¬ 
phologies on comparable angular scales. 


2.2.2 Large Millimeter Telescope AzTEC LI mm & Redshift 
Search Receiver 

The Large Millimeter Telescope Alfonso Serrano, hereafter LMT, 
is a 50m diameter millimeter-wavelength radio telescope on Vol- 
can Sierra Negra, Mexico, at an altitude of 4600 m. LMT ‘Early 
Science’ observations were conducted towards 9io9 with AzTEC 
(Wilson et al. 2008) and the Redshift Search Receiver (Erickson 
et al. 2007), operating in the 1mm and 3mm atmospheric windows 
respectively. During the ‘Early-Science’ programme only the inner 
32m diameter of the LMT primary surface was aligned and coupled 
to the scientific instruments. 

Continuum imaging of 9io9 was first conducted at 1.1mm with 
the AzTEC camera on 16th January 2014 UT, and then repeated on 
24th January 2014 UT due to a poor focus in the original observa¬ 
tions. The highest quality observations however were obtained dur¬ 
ing the early evening of November 2nd 2014 UT under improved 
weather conditions with a 225GHz zenith opacity of T 225 GHz ~ 



obs (GHz) 

Figure 3. The Redshift Search Receiver spectrum of 9io9 showing the 
full 38 GHz instantaneous bandwidth spectral coverage in the 3 mm atmo¬ 
spheric window (73-111 GHz). The detected molecular line at 97.312 GHz 
is ^^CO(3-2) from the lensed background source at a redshift 2 ; = 2.553. 

0.08. The measured beam-size at 1.1mm was 8.5" x 7.9" EWHM 
in azimuth and elevation respectively. The AzTEC observations 
were made in the standard Photometry Map Mode using a Lissajous 
scanning pattern, and reduced using the standard AzTEC data- 
analysis pipeline (Scott et al. 2008). Elux calibration was performed 
against Neptune, assuming a 1.1mm flux density of 18.38 Jy. A 
short four minute AzTEC integration detected the source 9io9 with 
a S/N=39 and a flux density of 95.5 zb 2.4 mJy, with no suggestion 
of extended emission. The accuracy of the absolute flux calibration 
is estimated to be ~10%. 

The LMT operates with an ultra-wideband spectrometer, the 
Redshift Search Receiver (RSR) that covers the full 3mm atmo¬ 
spheric window (73-111 GHz) at low spectral resolution (31,MHz 
channels, ~100km s“^ at 90 GHz). The RSR is a dual-beam dual¬ 
polarization spectrometer, with 38 GHz instantaneous bandwidth, 
designed to detect molecular gas, typically CO line transitions, 
from high-redshift galaxies. A 1 hour integration with the RSR 
was conducted towards 9io9 on 19th January 2014 UT with a mea¬ 
sured Tsys at 90 GHz of 85 K. The RSR data were reduced and 
calibrated using the LMT DREAMPY (Data REduction and Anal¬ 
ysis Methods in PYthon) package. A single strong spectral-line was 
detected at 97.31 GHz (Eigure 3). The combination of the photo¬ 
metric redshift of z ~ 2.6, determined from the (sub-)millimeter 
and EIR observations presented in this paper, with the RSR spec¬ 
trum eliminates alternative redshift solutions and any ambiguity in 
the measured spectroscopic redshift. We conclude that the detected 
line at 97.31 GHz is ^^CO(3-2) from the lensed background source 
at z = 2.553. Eull details of the LMT observations towards 9io9, 
together with a more complete analysis of the AzTEC and RSR 
data, are presented in Sanchez-Arguelles et al. (2015 in prepara¬ 
tion) and Harrington et al. (2015 in preparation). 

2.3 Near-infrared spectroscopy 

We observed 9io9 with the Infrared Camera and Spectrograph 
(IRCS, Kobayashi et al. 2000) instrument on the Subaru telescope 
on the night of 15 January 2014 UT as a back-up target for the 


© 0000 RAS, MNRAS 000, 000-000 




The Red Radio Ring 5 



1.4 1.6 1.8 2.0 2.2 2.4 

obs( m) 


Figure 4. Observed near-infrared HK spectrum of the primary arc of 9io9 
taken with IRCS on Subaru. The low resolution spectrum shows the clear 
rest-frame optical emission lines of the unresolved [O iii]A4959,A5007 
doublet, Ha+[N ii], and the [S ii] A6717,A6731 doublet, confirming the red- 
shift of the lensed galaxy at 2 ; = 2.553 (see §2.3). 
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Figure 5. Observed radio-to-sub-millimetre spectral energy distribution 
(SED) of 9io9, showing the integrated flux density of the lensed galaxy 
(we assume no contribution from the lensing galaxy at these frequencies). 
The SED is well-modelled by a combination of modified blackbody with 
T P:::^ 40 K, describing the cold dust thermal emission, and a simple power 
law at < 50 GHz, with spectral index a ^ 1.1, describing radio 
synchrotron emission. 


scheduled observing program^. The instrument was setup with 
the HK low-resolution grism eovering the wavelength range 1.4- 
2.45/im. We used the 0.9" wide slit of the ‘Reflective 3W’ slit- 
mask, resulting in a speetral resolving power of 77 95 at the 

central wavelength. The entranee slit was set at a position angle of 
155 degrees east of north to align it to the bright, western half of 
the red ring. We observed 9io9 using the nod-on-slit teehnique with 
slit positions separated by 7". We obtained 12 exposures for a to¬ 
tal of 1 hour on souree. Calibration data comprised of Argon lamp 
exposures for wavelength ealibration, and Halogen lamp exposures 
for flatfleld ealibration. 

The data were redueed with the IRAF twodspec package, and 
standard proeedures for flat-fielding and sky-subtraetion for the 
nod-on-slit teehnique were used. The final data produet is a two- 
dimensional speetrum of 9io9 with a regular sampling along the 
spatial and speetral direetions. To inerease the signal-to-noise ratio 
we binned the data in speetral direction by a factor of 6 x and in the 
spatial direction by a factor of 4 x. The total integrated speetrum of 
the western are of the red ring of 9io9 is shown in Figure 4. The 
redshifted nebular emission lines of Ha+[N ll] and (unresolved) 
[O lll]A4959,A5007 and [S ll]A6717,A6731 doublets are deteeted 
confirming the CO redshift of 2.553. 


3 ANALYSIS 

3.1 Spectral energy distribution 

In Figure 5 we plot the (integrated) observed radio-sub-millimetre 
speetral energy distribution (SED) of 9io9. We model the sub- 
millimetre-millimetre photometry with a single temperature (T) 
modified blackbody 


^ proposal ID S13B-002 


Su oc 


(1 - 

^hu/kT _ 


( 1 ) 


where the optical depth r{iy) = (y normalised to lOO/xm 
(rest) and /7 is the emissivity parameter. Sinee the redshift is aceu- 
rately known, we allow {3 and T to be free parameters. The best fit¬ 
ting parameters are [3 = 2.2zb0.4, T = (38zb4) K which provides a 
good model of the observed emission over 250-1 lOO/xm. The inte¬ 
grated (rest-frame) S-lOOO/xm luminosity is pLm = (1.3 zb 0.1) x 
10^^L©, where p is the lensing magnification. Note that this might 
underestimate the total infrared luminosity sinee we have not mod¬ 
elled mid-infrared emission from warmer dust components. We fit 
the radio portion of the SED with a power law Su oc with 
a = 1.11 zb 0.05. The eombined model SED is shown in Eigure 5. 


3.2 Lens modelling 

The system shows a primary lensing galaxy G1 and a seeondary 
satellite galaxy G2, probably at the same redshift (z = 0.206, see 
Eig. 1, 2). We ehoose an isothermal ellipsoid for the density pro¬ 
file of Gl. Sinee G2 is loeated elose to, and likely part of, the lens 
system, we model its density profile with an additional isothermal 
sphere. The parameters of the model are: positions and eritieal ra¬ 
dius (b) for Gl and G2, elliptieity (e) and position angle (PA) for 
Gl (^e), together with an external shear magnitude (7) and PA (0^). 
We use a Gaussian prior (e = 0.1 zb 0.1) on the elliptieity of Gl 
and a uniform prior on the position of G2 obtained from the optical 
light profile fitting of the galaxies. 

The lens mass model is optimized to simultaneously fit the ex¬ 
tended images from the FCg-band image (seeing 0 = (O.SzbO.l)") 
and the VLA 5 GHz data where the lensed are and the eounter im¬ 
age are deteeted with suffleient signal-to-noise. The absolute as¬ 
trometric calibration of the radio images is better than 100 mas, 
and we cannot improve the relative alignment of the near- 
infrared and VLA images. Furthermore, we acknowledge that 
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Figure 6. Mass modelling results based on joint fitting of the near-infrared imaging (colour map) and 5 GHz imaging (shown as contours at levels of 
0.0033 X (3,9,27) mJy beam” ^). (left) source plane model (centre left) the reconstructed image in the lens plane with the lensing galaxies G1 and G2 sub¬ 
tracted, (centre right) observed data linearly scaled between - 3 - 10 ( 7 , (right) residual (observed—model). The residual colour scale (near-infrared) is the same 
as the previous panel and the contours are at levels of 0.0033 x (-1,1) mJybeam”^. Images are orientated north up, east left. 


setting the source model to be coincident in the near-infrared 
and radio is an assumption that does not take into account the 
possibility that these have different source plane morphologies. 
However, it is necessary to make this assumption with the cur¬ 
rent data in order to reduce free parameters and degeneracies 
in the modelling. An annulus centered on the lens enclosing the arc 
and the counter image is used as a constraint when optimizing the 
lens model. We model the background source in the near-infrared 
with a Sersic profile (fiux, position, ellipticity and position angle, 
effective radius and Sersic index) and in the radio with a Gaussian 
profile. We set the sources at both wavelengths to be coincident as 
we expect them to be physically associated. 

To construct the mass model, we use gravlens (Keeton 2011) 
and a Markov Chain Monte Carlo (MCMC) simulation for the fit¬ 
ting procedure (using the Python module emcee, Foreman-Mackey 
et al. 2013) to optimize the parameters and explore degeneracies 
in the model. Initially, the position of G1 is held fixed during op¬ 
timization and after a reasonable model is found the position is 
allowed to be free. The source plane model, model lensed im¬ 
ages and residuals (for both wavelengths) are shown in Figure 6. 
The magnification in the near-infrared is /xnir ~ 11 and in the 
radio is /isGHz ~ 13, the unlensed magnitude of the source is 
Ks — 20.75 mag and S^ghz = 1.2 /xJybeam”^. The param¬ 
eters of the best fit lens model are summarized in Table 1. The 
critical radii from Table 1 can be taken with the lens and source 
redshifts to calculate the lens velocity dispersion for G1 and G2: 
cTGi = 316 zb 1 km s“^ and crG2 = 80 zb 4km s“^. 

Although the eMERLIN data are of somewhat lower signal- 
to-noise, their higher resolution allows at least preliminary investi¬ 
gation of the radio source structure. Using the galaxy mass model 
established above for the A^s-band and VLA imaging, we can at¬ 
tempt to reconstruct the eMERLIN data using a parametrised Gaus¬ 
sian fit to represent the source. The results are shown in Eigure 
7. In general, the source reconstruction is subject to considerable 
degeneracy because the main constraint provided by the lensing 
concerns the source distribution around the high-magnification re¬ 
gion corresponding to the caustic running across it. Nevertheless, 
the data indicate a double rather than a single Gaussian fit, if we 
believe the structure at the southern end of the arc which is diffi¬ 
cult to reproduce with a single component. The best-fitting model 
is subject to some degeneracies with the existing data, which can 
be investigated using a standard MCMC procedure. Specifically, 
the relative position of the centres of the components has an error 
approximately equal to their separation in the best-fitting model. 


Table 1. Best-fitting parameters describing the lens system comprised of 
components G1 and G2. Positions are given relative to the coordinates 
(02^09™41.268®, -1-00° 15'58.59"). Note that there is a degeneracy between 
the ellipticity of G1 and the shear magnitude and similarly, between the 
mass (or the critical radius) of G1 and G2. Uncertainties reflect the 16th and 
84th percentile range derived from the MCMC fitting of the mass model. 


Parameter 

G1 

G2 

X (mas) 

-94l?3 

37711^ 

y (mas) 

126lto 

233212^^ 

Critical radius (b, arcsec) 



Ellipticity (e) 

0.12 ± 0.02 

— 

PA (Oe, deg) 

84li 

— 

Shear magnitude ( 7 ) 

0.064lHo4 

— 

Shear PA (0^, deg) 

831^ 

— 


with the more compact component contributing most of the posi¬ 
tional error. The parameters describing the compact component, in¬ 
cluding its precise position, are uncertain apart from its size, which 
needs to be below approximately 30 mas, or 250 pc in the source 
plane. The more extended component is preferred to be elliptical 
{h/a — 0.271q iq with PA 15 zb 24 degrees), and together the radio 
structure is reasonably consistent with that of a small radio core and 
jet, although higher signal-to-noise data would be needed to make 
a more definitive statement, we defer this, together with a more 
systematic source reconstruction using other methods to a future 
work. 


4 DISCUSSION 

By simultaneously fitting the Kg-band and 5 GHz imaging we have 
derived a lens model that allows us to determine a magnification 
factor of approximately 10 x. The magnification-corrected infrared 
luminosity of 9io9 is Lir 10^^L©, ranking it amongst the most 
luminous galaxies in the Universe, in the ‘hyperluminous’ class. 
If the entire infrared luminosity budget is dominated by star forma¬ 
tion (i.e. re-processed ultraviolet emission from massive stars), then 
the implied star formation rate of 9io9 is M* 2500 M© yr“^ 
(Kennicutt & Evans (2012), assuming a Chabrier initial mass func¬ 
tion), a thousand times the star formation rate of the Milky Way. 
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Figure 7. (left) the reconstructed 1.4 GHz source model, using a two-Gaussian fit and the lens model described in §3.2. Note that the relative positions of 
the compact and extended components are not well constrained in the current model, (centre left) the reconstructed image in the lens plane, (centre right) 
eMERLIN observed data linearly scaled between -5-10(7, (right) residual (observed—model). The residual colour scale is the same as the previous panel. Note 
the scale bar for the source plane is 100 mas, corresponding to 820 pc in physical projection. With the exception of the first panel, images are centred on 
the position of the lensing galaxy G1 (Figure 1). The source plane image is slightly offset from centre for illustration purposes. Images are orientated 
north up, east left. 


Given the radio flux density and reconstructed radio morphology 
described above, it is likely that 9io9 contains an AGN that con¬ 
tributes non-negligibly to the energy output, and therefore M* 
should be taken as an upper limit, with the actual star formation 
rate scaling like 7 L 1 R, where 7 0.6 for starburst/AGN com¬ 

posite systems at similar epoch with comparable Lir, determined 
through mid-infrared diagnostics (Coppin et al. 2010). Even with a 
40% correction for AGN contribution to Lir, the implied star for¬ 
mation rate is still remarkably high, at over a thousand solar masses 
per year. 

In the absence of other diagnostics (such as the observation of 
broad lines, or emission line ratios such as [Olll]/H/3), the infrared- 
to-radio luminosity ratio can be used as a basic diagnostic for the 
presence of an active galactic nucleus, described by the parameter 

'Lir/3.75 X 


q = logi 


ii.4GHz/Wm-2Hz-i 


( 2 ) 


with the monochromatic radio luminosity L-corrected to the rest- 
frame 1.4GHz assuming the measured spectral slope. Star¬ 

forming galaxies have a rather tight distribution of q, with {q) ^ 
2.5 and scatter of approximately aq ^ 0.25 (e.g. Ivison et al. 
2010). The radio-far-infrared correlation is thought to arise from 
two linked processes associated with massive star formation: (a) 
cosmic-ray electrons from supernovae that emit synchrotron ra¬ 
diation as they are accelerated in magnetic Adds, leaking out of 
dust-enshrouded HII regions, and (b) infrared emission originating 
from dust-reprocessed ultraviolet photons emitted by young mas¬ 
sive stars. In galaxies hosting AGN, synchrotron emission associ¬ 
ated with black hole accretion can contribute to the radio flux, re¬ 
sulting in lower q. In this case, the magniflcation-corrected radio 
luminosity is Li. 4 ghz = (1.7 zb 0.1) x 10^^ WHz“^, implying 
^ = 1.8 zb 0.1 (here we have assumed the same magniflcation cor¬ 
rection for the radio and infrared emission, /x = 13). Thus, 9io9 is 
not consistent with a pure starburst, with signiflcant nuclear activity 
associated with the growth of the central black hole contributing to 
the radio emission. 

The current data allow us to postulate, with a reasonable de¬ 
gree of confldence, the basic nature of the radio source. There are 
some hints that the radio spectral index is flatter in the core com¬ 
ponent, a ^ 0.2 steepening to a 1 along the ring. This would 
be consistent with the classic picture of optically thick emission 
due to synchrotron self-absorption around a compact, dense nuclear 
source. The steepening spectrum could indicate aging of the rela¬ 


tivistic electron population (which lose their energy as L^) along 
a possible jet, but a more detailed analysis that takes into account 
the different visibilities sampled by eMERLIN and VLA, improved 
sampling of the radio spectrum and a more accurate lens model will 
be needed to investigate this. Considering the reconstructed source 
plane morphology, with compact core and potential jet/lobe struc¬ 
ture, the radio power of 9io9 places it in the Fanaroff-Riley II class 
of radio galaxies (FRII, Li . 4 GHz > W Hz“^). 


5 SUMMARY 

Through a citizen science project - Space Warps - we have discov¬ 
ered an extremely luminous galaxy at z = 2.5, close to the peak in 
stellar mass and black hole growth in the Universe, lensed by a fore¬ 
ground Luminous Red Galaxy, possibly part of a small group 
(given the large Einstein radius compared to typical single¬ 
galaxy lenses) at z = 0.2. The lensed source, dubbed 9io9, appears 
as a red Einstein ring, prominent in the near-infrared, but very faint 
in the optical (and thus would have been overlooked in an optical 
survey), characteristic of either a heavily dust-obscured or predom¬ 
inantly old stellar population. This lensed system would have been 
overlooked in an optical survey. Archival and subsequent follow¬ 
up observations revealed that 9io9 is extraordinarily bright in the 
sub-millimetre, with a flux density approaching 1 Jy at the peak of 
the thermal dust emission. We have resolved the lens at two radio 
frequencies, with sub-arcsecond eMERLIN and VLA mapping at 
1.4 GHz and 5 GHz respectively, revealing radio emission from the 
source tracing the near-infrared light. The strong radio and infrared 
emission indicates a galaxy undergoing an episode of rapid growth, 
either in an intense starburst episode, through signiflcant matter ac¬ 
cretion onto a central supermassive black hole, or both. 

A plausible evolutionary scenario for 9io9 is that the galaxy 
is destined to evolve into a massive (M* Mq) elliptical 

galaxy by 2 : = 0, perhaps located at the centre of a large group 
or cluster of galaxies. Seen when the Universe was less than 3 Gyr 
old, we are gaining an insight into the major growth phase of this 
massive galaxy, with gravitational lensing providing a unique win¬ 
dow onto spatial scales not possible for an equivalent non-lensed 
held source. The fact that high luminosity sources such as this are 
intrinsically rare even in the absence of lensing highlights the ben¬ 
efits that citizen science can offer in the mining of (ever growing) 
astronomical imaging archives, albeit at the cost of rather complex 
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selection functions for the objects found. We plan further follow-up 
of 9io9 to improve our understanding of the source. Particularly ex¬ 
citing is the possibility that we could detect and monitor variability 
in the source, which - if associated with the compact nuclear core 
- could be measured with cadences of days, weeks or months. 
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